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ABSTRACT 

Victoria-Regina isochrones for -0.4 < [a/Fe] < -1-0.4 and a wide range in [Fe/H], along with complementary 
zero-age horizontal branch (ZAHB) loci, have been applied to the color-magnitude diagram (CMD) of Carina. 

The color transformations that we have used have been “calibrated" so that isochrones provide excellent fits 
to the [(B-y)o, My]-diagrams of M3 and M92, when well supported estimates of the globular cluster (GC) 
reddenings and metallicities are assumed. The adopted distance moduli, for both the GCs and Carina, are based 
on our ZAHB models, which are able to reproduce the old HB component (as well as the luminosity of the HB 
clump) of the dwarf spheroidal galaxy quite well — even if it spans a range in [Fe/H] of ^ 1.5 dex, provided 
that [a/Fe] varies with [Fe/H] in approximately the way that has been derived spectroscopically. Ages derived 
here agree reasonably well with those found previously for the old and intermediate-age turnoff stars, as well 
as for the period of negligible star formation (SF) activity 6-10 Gyr ago). CMD simulations have been 
carried out for the faintest turnoff and subgiant stars. They indicate a clear preference for SF that lasted several 
Gyr instead of a short burst, with some indication that ages decrease with increasing [Fe/H]. In general, stellar 
models that assume spectroscopic metallicities provide satisfactory fits to the observations, including the thin 
giant branch of Carina, though higher oxygen abundances than those implied by the adopted values of [a/Fe] 
would have favorable consequences. 

Subject headings: galaxies: abundances — galaxies: dwarf—galaxies: individual (Carina) — galaxies: stellar 
content — globular clusters: individual (M 3, M 92) — stars: evolution 


1. INTRODUCTION 

The color-magnitude diagram (CMD) of the Carina dwarf 
spheroidal (dSph) galaxy has three prominent features that set 
it apart from those of similar systems. The first, and per¬ 
haps most conspicuous, characteristic is the remarkable di¬ 
chotomy of the horizontal branch (HB) into an intermediate- 
age, open-cluster-like red clump and an old, metal-poor, 
globular-cluster-like HB (see, e.g., Smecker-Hane et al. 1994, 
Hurley-Keller et al. 1998). In addition, deep photometric 
studies by Bono et al. (2010), Battaglia et al. (2012), and 
de Boer et al. (2014b), among others, have established that 
Carina has a very thin red-giant branch (RGB), despite the 
presence of stars that span a range in [Fe/H] from ~ -2.9 to 
^ -1.1 according to the spectroscopic surveys carried out by, 
e.g., Shetrone et al. (2003), Koch et al. (2008), Lemasle et 
al. (2012), and Venn et al. (2012). Finally, the former in¬ 
vestigations have shown that there is a clear gap in the dis¬ 
tribution of stars on the subgiant branch (SGB) that separate 
the oldest (> 10 Gyr) and intermediate-age (< 7 Gyr) stel¬ 
lar populations, and possibly on the main sequence (MS) at 
V ^ 22.5,B-V ^ 0.12 between stars with ages ^ 2 Gyr and 
those with higher ages. The evidence that Carina has under¬ 
gone at least two distinct, well-separated star formation events 
over its evolutionary history seems undeniable. 

Because Carina is so distant (> 105 kpc; e.g.. Vivas & Ma¬ 
teo 2013), high-resolution spectra have been obtained forrela- 
tively few of its stars, despite the very considerable effort that 
has been made by several groups (including those mentioned 
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above). As a consequence, the trends of [m/Fe] with [Fe/H] 
(for individual metals m) are not yet as well defined as they 
need to be in order for us to achieve a satisfactory understand¬ 
ing of Carina’s star formation history (SFH), and of its chem¬ 
ical evolution over time. Nevertheless, the available spectro¬ 
scopic results suggest that the lowest metallicity stars in Ca¬ 
rina have values of [a/Fe] Ri -1-0.4, as typically found in Milky 
Way halo stars (e.g., McWilliam 1997). However, in contrast 
with the Galaxy, where such enhancements are characteristic 
of most stars with [Fe/H] ^ -1.0, the transition to low values 
of [a/Fe] generally occurs at much lower iron abundances in 
dSph galaxies (see, e.g., Venn et al. 2004, Kirby et al. 2011). 
This is usually attributed to low star formation rates in the lat¬ 
ter (e.g., Gilmore & Wyse 1998), which would have the conse¬ 
quence that Type la supernovae (SN) could begin to contribute 
to their chemical evolution before the gas had been enriched 
to very high [Fe/H] values by exploding massive stars. 

In the case of the oldest Carina stars, the “knee" in the 
[a/Fe] vs. [Fe/H] relationship appears to occur at [Fe/H] 
^ -2.6 (de Boer et al. 2014a), where the plateau at lower 
metallicities changes into a roughly linear decrease of [a/Fe] 
to subsolar values at [Fe/H] ^-1.4. (Where this sequence 
terminates is not entirely clear due to the small numbers of 
observed stars and the considerable scatter in the derived 
abundances; also see de Boer et al. 2014b.) Regardless, be¬ 
cause the intermediate-age (lA) population includes stars that 
formed with [Fe/H] -1.8 and [a/Fe] -1-0.4, while the old¬ 
est population contains some stars at higher [Fe/H], it is clear 
that Carina has undergone inhomogeneous chemical evolution 
(Venn et al. 2012). This could plausibly be explained by gas 
infall (see the discussions by Lemasle et al. 2012 and de Boer 
et al. 2014b). (It is a reasonable expectation that each ma¬ 
jor SF event would produce an increase in [a/Fe] — but not 
necessarily to same value that is found in the preceding gen- 
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Figure 1. Pane/ (a): approximate mean [a/Fe] versus [Fe/H] relations for 
the oldest and intermediate-age stellar populations of Carina (the dashed and 
solid curves, respectively). The dashed curve, which reproduces the results 
given by de Boer et al. (2014a; their Fig. 3), clearly provides a good fit to the 
spectroscopic abundances compiled by Lemasle et al. (2012; see their Table 
10), assuming that [a/Fe] = [Mg/Fe] (see the text). The solid curve is sug¬ 
gested by a plot of the same data that was published by Venn et al. (2012, 
their Fig. 20). Only those stars in the Lemasle et al. sample that have sig¬ 
nificantly different iron and Mg, and hence a-element, abundances relative 
to the values along the solid curve have been plotted (the filled circles with 
attached vertical error bars) because we are mainly interested in the “Old" 
stars. Panel (b): schematic reproduction of the Caiina MDF that was derived 
by Starkenburg et al. (2010). The filled circles depict the variation in the 
number of stars with [Fe/H] that has been adopted in this study for the old¬ 
est Carina stars (those having chemical properties represented by the dashed 
curve) in the upper panel. 


eration — which subsequently falls with increasing [Fe/H] as 
the SF rate declines; see Gilmore & Wyse 1991.) 

Based on the Mg abundances reported by Lemasle et 
al. (2012, also see Venn et al. 2012), the Mg/Fe number abun¬ 
dance ratios in stars belonging to the lA population appear to 
decrease from [Mg/Fe] « -1-0.4 at [Fe/H] -1.8 to ~ 0.0 at 
[Fe/H] ~ -1.2 (but with a large scatter, mostly to very low val¬ 
ues of [Mg/Fe]). The same studies show that calcium abun¬ 
dances are considerably more uniform, suggesting no more 
than a shallow slope in the [Ca/Fe] vs. [Fe/H] relation over 
the same range in iron content. Because magnesium has a 
far bigger impact on stellar models due to its much greater 
absolute abundance (see VandenBerg et al. 2012, hereafter 
V12), Mg is arguably a better tracer of [a/Fe] than Ca. For 
this reason and the above considerations, we have adopted 
the relations between [a/Fe] and [Fe/H] for the oldest and 
lA populations of Carina that are shown in Figure la. The 
abundances for individual stars have been plotted only if they 
have [Fe/H] < -1.8 or they have [Mg/Fe] values that differ by 
more than 0.25 dex, at their [Fe/H] values, from those indi¬ 
cated by the solid curve. These data suggest that the dashed 
locus probably does extend to [Fe/H] '^-1.2 (but see de Boer 
et al. 2014b). Furthermore, the observed large scatter in the 
Mg, and hence a-element, abundances at [Fe/H] ^ -1.8 is a 
natural outcome if stars with these iron abundances lie along 
two, well-separated [a/Fe] versus [Fe/H] relations. 

The Carina metallicity distribution function (MDF) that 
was derived by Starkenburg et al. (2010), on the basis of their 
improved calibration of Ca II triplet absorption line strengths 


in terms of [Fe/H], is given in Fig. lb. Most of the stars in Ca¬ 
rina have [Fe/H] values between ^ -2.8 and -1.3, with just a 
few stars at both lower and higher metallicities. Local peaks 
in the number count distributions occur at [Fe/H] Ri -2.2 and 
-1.65. Although based on a much smaller sample, the MDF 
derived more recently by Lemasle et al. (2012) agrees quite 
well with the Starkenburg et al. findings. The main difference 
between them is at [Fe/H] ~ -2.0, where significantly fewer 
stars, relative to the peaks on either side, are present in the 
Lemasle et al. sample. 

Insofar as ages are concerned, the narrow RGB of Carina 
together with the lack of any obvious metallicity dependence 
of the colors of red giants at a common luminosity implies 
that ages and metallicities are anti-correlated (see, e.g., Koch 
et al. 2006)0 Although there have been some attempts to de¬ 
termine the ages of bright giants using isochrones for close 
to the spectroscopically derived [Fe/H] values (e.g., Lemasle 
et al. 2012), “the small offset in color with age ... combined 
with the sensitivity to the largely unknown elemental abun¬ 
dance mix and He content limits the usefulness of this ex¬ 
ercise" (Koch et al. 2006). Any dispersion in the observed 
abundances of the a-elements at the same [Fe/H] will have 
important consequences for the effective temperatures (and 
presumably the colors) of red giants, at a fixed age (as shown 
by V12). Ages are much more reliably determined for turnoff 
(TO) and SCjB stars since, at low metallicities, TO luminosity 
vs. age relations depend almost entirely on the absolute abun¬ 
dance of oxygen (and of C and N); the other metals mainly 
affect the predicted Teff scale in a relatively minor way (see 
V12). Unfortunately, as noted by Monelli et al. (2014), very 
little is known about the abundances of the CNO elements, 
and their variations with [Fe/H], in this (or any other) dSph 
galaxy. 

When this project began (in mid-2014), it was our inten¬ 
tion to use the latest Victoria-Regina isochrones (VandenBerg 
et al. 2014a; hereafter V14a) for -0.4 < [a/Fe] < -1-0.4 (and 
a wide range in [Fe/H]), together with new, complementary 
color transformations (Casagrande & VandenBerg 2014, here¬ 
after CV14), to study the implications of the very narrow RGB 
of Carina, and to try to shed some light on a possible con¬ 
flict between photometrically inferred and spectroscopically 
derived metallicities. A few years earlier, Bono et al. (2010; 
also see Lianou et al. 2011) had concluded, mainly from the 
superposition of Galactic star cluster fiducial sequences onto 
the Carina CMD and theoretical simulations of the latter, that 
the spreads in [Fe/H] encompassed by both the oldest and the 
intermediate-age (lA) populations in the dSph galaxy appear 
to be considerably less than those derived from spectroscopic 
data (also see Fabrizio et al. 2012). 

However, there is a third motivation for the present under¬ 
taking. When interpreting the photometry of complex stellar 
populations (as in the case of Carina), it is important to check 
that the isochrones which are used to determine their proper¬ 
ties (e.g., ages, SFH) faithfully reproduce the CMDs of sim- 
ple(r) systems, like globular clusters (GCs). If they fail to do 
so, they cannot be expected to overlie stars in galactic CMDs 
that have the same age and chemical abundances, but different 
evolutionary states (due to differences in mass). In fact, it is a 

Only recently has it been possible to separate the most metal-poor gi¬ 
ants from those having intermediate metallicities using measurements of the 
= (U index (Monelli et al. 2014). However, as discussed in 
the latter study, it has proven to be difficult to explain such obseiwations using 
current isochrones, perhaps because the assumed heavy-element mixtures are 
not representative of the actual stellar abundances. 
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good idea to “calibrate" stellar models using GCs and/or open 
clusters before applying them to complicated systems — as 
in the study of Andromeda’s SFH by Brown et al. (2006) and 
their recent analyses of the CMDs for several ultra-faint dwarf 
(UFD) galaxies (Brown et al. 2014, 2012). This has gener¬ 
ally not been done in studies of Carina, but we have opted to 
do so here. [Note that the ages of most star clusters are not 
subject to this concern because they can be derived by match¬ 
ing isochrones to just a small part of the CMD containing the 
turnoff luminosity (via, e.g., the procedure described by Van- 
denBerg et al. 2013; hereafter V13). In the case of galaxies, 
ages are effectively determined on a star-by-star basis, and the 
isochrones which are used should provide good fits to photo¬ 
metric data in all of the evolutionary stages that are studied in 
order to obtain the most reliable results.] 

The next section describes the observations that are em¬ 
ployed in this investigation, both of Carina and of the GCs 
M3 and M92, which are used to calibrate the stellar mod¬ 
els (see § 3) that have been applied to the CMD of the dSph 
(as reported in § 4). Because the most recent Victoria-Regina 
isochrones (V14a) are restricted to ages ^ 5 Gyr, the main 
focus of this project is on the faintest TO stars in Carina, 
which we have attempted to explain through the superposition 
of isochrones onto their observed CMD locations, as well as 
numerical simulations. The latter enable us to probe the SFH 
of Carina. A short summary of our findings is given in § 5. 

2. PHOTOMETRY OF CARINA, M 92, AND M 3 

The data for Carina discussed here are the same as those 
that were the subject of the analysis by Bono et al. (2010). 
They resulted from the analysis of more than 4,000 CCD im¬ 
ages obtained during 27 observing runs between 1992 De¬ 
cember and 2008 September. The observational indices from 
each observing run were transformed to the standard John¬ 
son UBV, Kron-Cousins RI photometric magnitude system of 
Landolt (1992; supplemented by additional measurements in 
the UBV filters by Landolt 1973) following the procedures 
described by, e.g.. Stetson (2000, 2005). Then the average re¬ 
sults from all observing runs were used to define a network of 
local photometric standard stars in the Carina field, to which 
observations of all other stars in the field were referenced. 

Exactly the same procedures were employed to produce 
the photometric results for our comparision clusters, M 3 (= 
NGC5272) and M92 (= NGC6341). The results for these 
two clusters come from 21 observing runs (1984 June-2002 
March) for M 3 and 40 observing runs (1984 June-2002 May) 
for M 92. 

In the final analysis, these results are tied to 324 and 336 
Landolt stars in the B and V filters, and 231 Landolt stars in 
the / filter. If we consider the V filter as representative, the 
median number of observations per star by Landolt is 16; the 
median number of observations per star by us is 69. We find 
that the r.m.s. dispersion of the difference between Landolt’s 
and our average photometric results is 0.012 mag per star in 
B, 0.009 mag in V, and 0.011 mag in I. We believe that this 
represents an irreducible limit to our ability to transform in¬ 
dices from one typical set of broad-band filters to another, 
given the inherent differences in the detailed spectral energy 
distributions of stars of a given perceived color. The net dif¬ 
ference between Landolt’s overall photometric system and our 
attempt to reproduce that system should be of order 0.01 mag 
divided by the square root of a number of order a few hundred, 
at least within the range of temperature, gravity, and chemical 
abundance that is spanned by Landolt’s (mostly Population I) 


photometric standards. 

It is true that our photometric study of Carina requires that 
we apply transformations derived from field standard stars of 
mostly solar metallicity to program stars that are apprecia¬ 
bly metal-poor. To the extent that different filter sets dif¬ 
ferently sample the stellar spectral-energy distributions, this 
could in principle introduce a systematic difference between 
the results of any one of our filter sets and Landolt’s filter set. 
Given that our 27 observing runs employed at least six differ¬ 
ent combinations of detector and filter set, this should miti¬ 
gate any net systematic difference between the average pho¬ 
tometric indices obtained by us and those that Landolt would 
have obtained by observing Carina with his equipment. In 
any case, if there does remain a metallity-dependent differ¬ 
ence between what we did measure and what Landolt would 
have measured, it should affect Carina, M3, and M92 sim¬ 
ilarly, and should be negligible in differential comparisons. 
We therefore believe that comparisons between our observa¬ 
tions and our theoretical predictions will be more hampered 
by imperfect approximations in the theoretical modelling, by 
uncertain transformations between the theoretical and obser¬ 
vational indices, and perhaps by imperfect knowledge of the 
reddening toward Carina, than by random or systematic errors 
in the photometry. 

To select our samples of well-observed, probable members 
of each system, we first estimated the photometric barycenter 
of the target and then plotted the relationship between appar¬ 
ent magnitude and photometric uncertainty as a function of 
the radial distance of the star from the center of its host sys¬ 
tem. In particular, we plotted V-band magnitude versus radius 
for stars with 0.01 < a{B-I) < 0.02mag, 0.04 < a{B-I) < 
0.06mag, and 0.08 < aiB-I) < O.lOmag. We then chose an 
annulus where all of these radial plots were essentially flat. 
Toward smaller radii the errors increase (or, more correctly 
stated, the apparent magnitude typical of a given photometric 
error becomes brighter) due to crowding; toward larger radii, 
the errors increase because those parts of the field typically 
appear in fewer CCD images. As a further limit on the outer 
radius of the annulus we plotted observed color versus radius 
for stars within a small apparent magnitude range around the 
cluster turnoff. We made sure that stars near the turnoff color 
dominated the plot out to our adopted outer limit. Lor M 3 
we adopted an annulus 120 < r < 360 arcseconds, while for 
M92 we adopted 120 < r < 600 arcseconds. Lor Carina we 
adopted the entire disk r < 600arcsec because we detected no 
inward increase in photometric scatter at any radius; appar¬ 
ently, crowding is not a limiting factor in our data for stars of 
any brightness within the magnitude range of interest at any 
radius all the way in to the center of the galaxy. 

3. STELLAR EVOLUTIONARY MODELS 

Most of the isochrones considered here were derived from 
the grids of evolutionary tracks that have been made pub¬ 
licly available by V14a, using their interpolation software, 
and converted to the observed planes via the transformations 
produced by CV14. Both the stellar models and the color- 
Teff relations can be interpolated to arbitrary a-element abun¬ 
dances within the range -0.4 < [a/Le] < -1-0.4 at nearly all 
of the [Le/H] values of interest. (Lor this study, we have as¬ 
sumed Y = 0.25 for the initial mass-fraction abundance of he¬ 
lium, though it is possible to generate isochrones from the 
V14a grids for any value of Y between 0.25 and 0.33.) 

We have also made use of a few supplementary grids for 
[Le/H] < -2.4 that will be the subject of a forthcoming pa- 
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per on extremely metal-deficient stars (those having values 
of [Fe/H] as low as -4.0). Some of these computations, 
which allow for variations in [O/Fe] from -tO.4 to -tl.2 at 
fixed [Fe/H] values, assuming [m/Fe] = -tO.4 for the other 
a-elements, have been employed by Brown et al. (2014) to 
interpret their Hubble Space Telescope (HST) photometry of 
several UFD galaxies. In addition, we are in the process of 
generating fully consistent zero-age horizontal branches (ZA- 
HBs) for all of the aforementioned grids and se veral of them 
appear in the plots that are presented in Section |4T| (One of 
the V14a grids and three ZAHBs have been extended to suf¬ 
ficiently high masses to permit comparisons of stellar models 
with the observations of the intermediate-age TO stars in Ca¬ 
rina and with the location of its red HB clump.) 

3.1. Calibration of the Isochrones 

As discussed at some length by V13 (also see V14a), there 
are so many uncertain factors that play a role in intercompar¬ 
isons of predicted and observed CMDs (such as the distance 
and [Fe/H] scales, heavy-element mixtures, color transforma¬ 
tions, convection theory, atmospheric boundary conditions, 
the role of diffusive and extra mixing processes) that “perfect" 
agreement between observations and theory is not a realistic 
expectation. It is, in fact, remarkable that current isochrones 
fare as well as they do, as the discrepancies that are found in, 
for instance, fits of current Victoria-Regina stellar models to 
GC photometry (see V14a, their Fig. 14) generally involve no 
more than relatively small zero-point and systematic offsets 
when well-supported estimates of the cluster properties are 
assumed. These offsets are typically at the level of ^ 0.025 
mag in the case of V-Zc colors (see CV14) and the very sim¬ 
ilar colors which can be derived from the HST FbQbW and 
F814VT bandpasses (VI3, V14a). 

For the most part, isochrones seem to perform nearly as 
well on the (B-V, V)-diagram as on the (V-Ic, y)-plane (see 
the examples given by CV14). However, as discovered in the 
course of this study, B-V colors appear to become increas¬ 
ingly problematic at the lowest [Fe/H] values. Whereas the 
isochrones provide rather good fits to our BV data for M 3 
(comparable to those presented for M5 by CV14) when rea¬ 
sonable assumptions are made concerning its basic proper¬ 
ties (i.e., distance, reddening, and metallicity), the agreement 
is noticeably less satisfactory in the case of M92. While 
the discrepancies are relatively minor above the turnoff, the 
models (for [Fe/H] = -2.40; Kraft & Ivans 2008; Carretta et 
al. 2009a) predict an appreciably steeper MS slope than the 
observed one. It seems doubtful that this is due to problems 
with the model Teff scale because the isochrones do not have 
similar difficulties when confronting M 92 data derived from 
longer wavelength filters (e.g., V -Ic, or the HST equivalent; 
see VI3, V14a). It is even less likely that this difficulty can be 
attributed to observational errors because our CMD for M 92 
is based on a very large number of independently calibrated 
data sets (see §|2]l. 

However, a resolution of this issue is not needed for this 
project. Since the BV observations of Carina, M3, and M92 
have all been calibrated to the same photometric system, we 
can force the stellar models to be on very close to the same 
system if we simply correct the colors along the isochrones so 
that they reproduce the CMDs of the two GCs. Even though 
this does tie our analysis to particular distance and [Fe/H] 
scales for globular clusters, the resultant models should yield 
a more accurate interpretation of the Carina photometry than 
those without such adjustments. One must just be careful to 


carry out this calibration of the synthetic colors in such a way 
that the latter vary smoothly with [Fe/H] (as well as with Teff 
and log g) in order that they can be applied to isochrones for 
the ranges in age and metallicity found in Carina. 

We believe that the best way to achieve this is to set up 
tables of 5{B-V) values to complement the tabulated trans¬ 
formations to B-V given by CV14. As these results will 
be used primarily in the fitting of isochrones to the faintest 
TO of Carina, the color correction tables can be restricted 
to [Fe/H] =-3.0(0.5)-1.5, and at each [Fe/H] value, to 
logg = 2.5(0.5)5.0 and T^« = 5000(250)6750 K, where the 
numbers in parentheses give the increments to the three vari¬ 
ables. The color-Teff relations provided by CV14 consider 
exactly the same values of [Fe/H], logg, and T^s among their 
coverage of a greatly expanded parameter space (and ~ 40 
different color indices). They also provide tables for [a/Fe] 
= -0.4(0.4)-t0.4, but the majority of the Galactic GCs appear 
to have [a/Fe] ss -1-0.4 (see, e.g., Carretta et al. 2009b); con¬ 
sequently, the S(B-V) adjustments that are derived from a 
consideration of the M 3 and M 92 CMDs implicitly assume 
this value of [a/Fe]. (In fact, we have not attempted to de¬ 
termine how 6{B - V) depends on the abundances of the a- 
elements. As shown below, these corrections are, in any case, 
quite small for the post-turnoff portions of Victoria-Regina 
isochrones, and since they are used in conjunction with the 
CV14 transformations, which do predict the effects on B-V 
of differences in [a/Fe], this neglect should not be a concern. 
In other words, we have assumed that, at a fixed value of 
[Fe/H], the dependence of 6{B-V) on logg and Teff is not 
a function of [a/Fe].) 

M 3 and M 92 are particularly suitable GCs for this exer¬ 
cise because their [Fe/H] values are close to the grid values of 
-1.5 and -2.5, respectively. Hence, whatever adjustments in 
the B-V color are needed to obtain good fits of isochrones to 
the CMDs of these systems will be very similar to the values 
that should appear in the interpolation tables for these metal- 
licities. Because B-V is much more dependent on Teff than 
on log g (especially for MS, SGB, and lower RGB stars; see, 
e.g., VandenBerg & Clem 2003, their Fig. 3), the temperature 
shift, at a fixed gravity, between the MS portions of isochrones 
for [Fe/H] = -2.0 and -2.5, as a fraction of that predicted for 
[Fe/H] =-1.5 and -2.5 isochrones of the same age, should be 
a reasonable approximation to the corresponding variation in 
S(B-V). Based on our stellar models, this fraction is close 
to 0.30, which we have adopted in the creation of a table 
of color corrections for [Fe/H] = -2.0 (and for -3.0, using 
a similar rationale). Indeed, with relatively few iterations, it 
was quite easy to produce a table of 6{B - V) values that are 
smooth functions of [Fe/H], logg, and Teff and which enable 
isochrones to satisfy the M 3 and M 92 contraints very well. 
Importantly, these corrections are small for the age-sensitive 
parts of isochrones — easily within the uncertainties associ¬ 
ated with the stellar models, color-Tiff relations, and photo¬ 
metric zero-points. 

This can be seen in Figure 2, which compares the transfor¬ 
mations to B-V for [Fe/H] = -2.0 (and [a/Fe] = -tO.4) from 
CV14 (filled circles) with those obtained if the latter are cor¬ 
rected by our adopted values of 6(B-V) (open circles). A 
12 Gyr isochrone for the same chemical abundances has also 
been plotted to permit visual estimates of the variations in 
the color adjustments that have been applied as a function of 
evolutionary state. Above the turnoff, they amount to < 0.02 
mag, rising to ~ 0.04 mag only at log g ^ 4.7. Lines of con¬ 
stant Teff run in a nearly vertical direction; e.g., the bluest 
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Figure 2. Plot of the (5-y)-reff-log g relations for [Fe/H] = -2.0, log g = 
2.5(0.5)5.0 and, at constant gravity, Feff = 6750(-250)5000 K, in the direction 
from left to right. The filled and open circles give, in turn, the transformations 
reported by CV14 and those derived here from fits of isochrones to the M3 
and M92 CMDs, and other considerations (see the text). The S{B—V) values 
refemed to in the text correspond to the differences bet\veen the open and 
filled circles at common values of log g and Teff. The solid curve represents a 
12 Gyr isochrone for [Fe/H] =-2.0, [ci/Fe] =+0.4, and Y = 0.25 (from V14a). 
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Figure 3. Fits of isochrones for the indicated ages and chemical abundances 
to our CMDs for M92 and M3 on the assumption of the apparent dis¬ 
tance moduli (based on ZAHB models) and reddenings (from Schlafly & 
Finkbeiner 2011) that are specified in the upper left-hand comers of the re¬ 
spective panels. The predicted B-V colors were derived from the CV14 
transformations, but with small adjustments to the latter that were obtained 
by interpolation in tables of 5{B-V) values for [Fe/H] = —3.0,-2.5,-2.0, and 
-1.5 and ranges of log g and Tgff relevant to the faintest TO stars in Carina. 
(The differences between the adopted and CV14 color-Teff relations for stars 
having [Fe/H] = —2.0 are illustrated in the previous figure.) 


points at the grid values of log g from 2.5 to 5.0 give the pre¬ 
dicted colors at T^ff = 6750 K, whereas the reddest ones indi¬ 
cate the colors for the coolest T^ff considered (5000 K). At a 
fixed gravity, the points at intermediate B-V values illustrate 
how this color index varies when the temperature decreases in 
successive steps of 250 K, in the direction from left to right. 

Figure 3 demonstrates that our calibration of the CV14 


Figure 4. Plot of isochrones for the indicated age and chemical abundances if 
the modified CV14 transformations described in this paper ai‘e used to trans¬ 
pose them to the (B - V, Mv)-diagram. 


transformations does, indeed, lead to excellent fits of 
isochrones to the CMDs of M3 and M92 (by design, of 
course), when the reddenings given by Schlafly & Finkbeiner 
(2011) and apparent distance moduli based on ZAHB mod¬ 
els for the indicated chemical abundances are assumed. (Dis¬ 
tances derived in this way are in very good agreement with 
those based on the RR Lyrae standard candle; see V13 and 
VandenBerg et al. 2014b, hereafter V14b.) It is worth point¬ 
ing out that the same ages would be (and are) found using the 
CV14 color-Feff relations because the predicted TO luminosi¬ 
ties have not been altered by the adopted d(B-V) offsets. The 
main effect of the latter, besides removing small discrepancies 
along the lower RGB, is to force the predicted and observed 
TO colors (and lower main sequences) to agree. This concur¬ 
rence is required, in fact, if one is to obtain the best estimate 
of the turnoff age of a GC from fits of isochrones to just the 
small region of its CMD in the vicinity of the TO (see the 
discussion of this issue by V13, their § 3.1)0 In any case, the 
colors along Victoria-Regina isochrones should he adjusted in 
the way that we have described if the ZAHB-based distance 
scale is accurate and the actual cluster [Fe/H] values are close 
to the values that we have adopted (see Fig. 3). 

To further illustrate the performance of the calibrated trans¬ 
formations to B-V, we have plotted, in Figure 4, isochrones 
for seven [Fe/H] values (as indicated), each assuming Y = 
0.25, [a/Fe] = -1-0.4, and an age of 12 Gyr. The differences be¬ 
tween them, at a fixed value of My, vary smoothly and mono- 
tonically, and there are no obvious irregularities as a function 
of absolute magnitude or color. The main effect of the S{B-V) 

^ The main reason why V13 obtained ages for M 3 and M 92 that are 0.5 
Gyr younger than those derived here — despite their assumption of slightly 
smaller distance moduli by 0.03 mag, which would normally imply increased 
ages by ^ 0.3 Gyr — is that the model grids used in their study assumed 
higher oxygen abundances by 0.24 dex (in an absolute sense), due to differ¬ 
ences in the respective solar heavy-element mixtures and the adopted values 
of [0/Fe]. The higher the absolute oxygen abundance, the lower the age at a 
given TO luminosity (see V12). This example serves to illustrate the impor¬ 
tance of taking chemical composition and distance modulus differences into 
account when comparing age determinations for the same star cluster that are 
reported in different papers. 
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values that we have generated is to displace the lower-MS por¬ 
tions of the isochrones to redder colors, especially at [Fe/H] 
;$ -1.8, by amounts that vary directly with My- As a result, 
the isochrones display a considerably reduced sensitivity to 
[Fe/H] over the range 4 < My ^ 6, as required by the pho¬ 
tometry of M 3 and M 92 (if our assumed properties for them 
are accurate), compared with their behavior when the CV14 
transformations are used without any adjustments. 

Figs. 2, 3, and 4 thus give us considerable confidence that 
the isochrones used in this investigation to model the oldest 
MS, SGB, and lower RGB stars in Carina are about as well 
constrained to empirical data as we are able to make them. 
In the case of giant stars (those with logg < 2.5), the CV14 
transformations to B-V appear to require no more than a 
-0.01 mag adjustment in order for the same isochrones that 
have been plotted in Fi g. 3 to reproduce the entire RGBs 
of M 3 and M 92 (see § 14.31 which presents an analysis of 
the giant-branch populations of the dSph). Finally, we note 
that the interpolated 6(B — V) values for Feff = 6750 K are 
used when the predicted temperatures exceed this limit (as 
in isochrones that are relevant to lA stars), and the tabulated 
color adjustments for [Fe/H] = -1.50 are adopted for slightly 
more metal-rich stars. The latter is a reasonable assumption 
given that the CV14 transformations appear to provide com¬ 
parable fits of isochrones to the CMD of M5 (see CV14), 
which has [Fe/H ~ -1.3 (e.g., Carretta et al. 2009a), and of 
the ~ 0.2 dex more metal-deficient cluster M 3, as we found 
during the course of this work. Moreover, there appear to be 
very few stars in the faintest TO of Carina with [Fe/H] >-1.5 
(see Fig. la); consequently, they will have no more than a very 
minor influence on our findings. 

4. INTERPRETATION OE THE CARINA CMD 

To try to constrain our understanding of Carina, we will 
hrst intercompare its CMD with that of M92. In addition to 
highlighting some of the similarities and differences between 
these systems, this will provide a useful check of the reliabil¬ 
ity of the adopted reddenings and distance moduli in a relative 
sense. This will be followed by hts of isochrones that utilize 
the empirically calibrated color transformations just described 
to the CMD of the galaxy, to explore the implications of dif¬ 
ferent assumptions concerning the ages and chemical compo¬ 
sitions of its stellar populations. Eully consistent ZAHB loci 
will also be considered. Whereas the analyses presented in 
the first subsection will reveal some of the possible interpre¬ 
tations of Carina’s CMD, with particular emphasis on the old¬ 
est stars, the next one will endeavor to assess the viability of 
these possibilities through numerical simulations. What does 
the observed MDE and the luminosity width of the SGB asso¬ 
ciated with the faintest TO tell us about the SEH? The third, 
and final, subsection will examine the thin RGB of Carina to 
ascertain the extent to which it is the outcome of a fortuitous 
alignment of the ages and metallicities of the stars that popu¬ 
late this feature. 

4.1. Turnoff and HB Stars in Carina 

In the upper panel of Eigure 5, the CMD of M92 has 
been transposed to the [(B-y)o,My]-diagram assuming (as 
adopted previously; see Eig. 3) E(B-V) = 0.019 and (m- 
M)y = 14.75. This apparent distance modulus is obtained if 
a ZAHB for Y = 0.25 and [Ee/H] = -2.40 (the solid curve) 
is htted to the lower bound of the distribution of the clus¬ 
ter HB stars (as shown). To obtain this result, the synthetic 


B-V colors along the ZAHB were adjusted to the red by 
0.016 mag, which is obtained from our S(B-V) tables if 
[Ee/H] = -2.40 and Teff > 6750 K. Interestingly, V14b (see 
their Eig. 5) found that the same ZAHB provided an equally 
good fit to the HB of M 92 on the [(V-/c)o, Mi/]-diagram, but 
without having to apply any adjustment of the synthetic colors 
(from CV14). Essentially the same thing was found by V13 
(see their Eig. 11) in a study of the nearly equivalent CMD 
of M 92 that they constructed from the HST F6Q6W, F’814W 
observations that were obtained by Sarajedini et al. (2007) 
— though the stellar models employed in the 2013 investi¬ 
gation assumed a somewhat different metals mixture. [While 
the main cause of such minor inconsistencies is difficult to de¬ 
termine, it should not be a surprise that they exist given that, 
e.g., on-going improvements to model atmospheres and the 
line lists used in synthetic spectra are bound to affect the mag¬ 
nitudes derived from some filter bandpasses more than others. 
In particular, B is affected by line blanketing more than V and 
/c; consequently, B-V will be a more challenging color index 
to model than V-ff-] 

If the reddening in the direction of Carina is E{B-V) = 
0.054 (Schlafly & Einkbeiner 2011), the adoption of (m- 
M)y = 20.28 results in the superposition of the Carina and 
M92 CMDs shown in the top panel. In this case, the cluster 
SGB is coincident with the location of the densest concentra¬ 
tion of Carina subgiants that are associated with the faintest 
TO, implying a similar age if the latter have the same metal- 
licity. However, judging from the very different luminosity 
widths of their respective SGBs, the oldest stellar population 
in Carina must span a much wider range in age or metallicity, 
or both, than the stars residing in M92. (These indications 
from photometry have, of course, been confirmed by several 
spectroscopic studies; see §[T]) The assumed reddenings and 
distance moduli also imply that the RGB and HB of M 92 lie 
along the blue edges of the galactic distribution of giants and 
horizontal branch stars, respectively, which is consistent with 
M92 being more metal-poor than the majority of the stars 
in Carina. Encouragingly, the adopted apparent modulus im¬ 
plies a true distance modulus of (m-M)o « 20.11, which is 
in very good agreement with recent determinations based on 
near-IR magnitudes of the RGB tip (20.08-20.12, Pietrzyhski 
et al. 2009), on the RR Lyrae variables, (20.09, Coppola et 
al. 2013), and on the dwarf Cepheids in Carina (20.17, Vivas 
& Mateo 2013). 

The location of the ZAHB for [Ee/H] = -2.4 relative to 
the blue HB stars in Carina is more easily seen in the bot¬ 
tom panel, along with the same isochrone (for 13 Gyr) that 
matches the TO and SGB of M92. This isochrone lies below 
most, but not all, of the faintest subgiants in the dSph: a bet¬ 
ter match to the least luminous SGB stars is clearly provided 
by a 13 Gyr isochrone for [Ee/H] = -2.0. However, this is no 
more than suggestive, given that the assumed value of [a/Ee] 
in these models is too high (recall Eig. la), and because a 
similar ht could be obtained using an even more metal-rich 
isochrone for a somewhat younger age. Still, these computa¬ 
tions give the impression that Carina contains stars at least as 
old as those found in the most metal-poor GCs (^13 Gyr). 
In addition, they illustrate that the luminosity width of Ca¬ 
rina’s SGB is consistent with a range in age ^ 4 Gyr, at a 
fixed [Ee/H] (note the separation between the 9 and 13 Gyr 
isochrones for [Ee/H] = -2.0), or a variation in metallicity, at 
a given age, that is much larger than 0.4 dex (note that the 
13 Gyr isochrones for [Ee/H] = -2.0 and -2.4 enclose only a 
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Figure 5. Upper panel: Superposition of the M 92 CMD onto that of Caiina 
assuming the indicated values of the reddening and apparent distance moduli. 
The solid curve represents a ZAHB for [Fe/H] = -2.40, [ct/Fe] = +0.40, and 
Y = 0.25. Lower panel: Overlay of isochrones for the indicated chemical 
abundances and ages, along with fully consistent ZAHB loci, onto the CMD 
of Carina. Only the red end of the ZAHB is shown for [Fe/H] = —1.8 and 
[a/Fe] = +0.4 (long dashes), as well as the one for [Fe/H] = —1.4 and [a/Fe] 
= 0.0 (the curve with a small filled circle attached to the bright end). 

small fraction of the observed subgiants). 

The other isochrones which have been plotted indicate that 
little or no star formation occurred in a 2-3 Gyr period prior 
to the major SF burst that happened about 6 Gyr ago, and that 
there may have been a lull in the SF activity ^2.5 Gyr ago as 
there is some (albeit marginal) indication of enhanced num¬ 
bers of MS stars at both younger and somewhat older ages. It 
is worth mentioning that, if isochrones for Y = 0.25, [Fe/H] 
= -1.4, and [a/Fe] =-1-0.13 (parameters that are also con¬ 
sistent with the [a/Fe] vs. [Fe/H] relation shown in Fig. la) 
were compared with the Carina CMD, we would obtain an 
age of 5.9 Gyr for the oldest of the intermediate-age stars 
(as compared with 6.3 Gyr if they have [Fe/H] = -1.8 and 
[a/Fe] = -1-0.4; see Fig. 5). The metallicity dependence of this 
age estimate is thus quite small — because the change in the 
absolute oxygen abundance between these two cases is only 



(B-V)„ 

Figure 6. Similar to the bottom panel in the previous figure, except that 
isochrones and ZAHBs for values of [Fe/H] and [a/Fe] that are more consis¬ 
tent with spectroscopic determinations (see Fig. la) ai‘e compared with the 
Carina CMD. 

0.13 dex. As discussed by V12, isochrones for the same oxy¬ 
gen abundance (i.e., [O/H]), but different [Fe/H], will predict 
nearly identical TO luminosity vs. age relations (especially at 
low metallicities, where the contributions of the metals to the 
opacities in the interiors of stars are small). 

Aside from a 1-2 Gyr reduction in the absolute ages, these 
results are quite similar to those reported by, e.g., Hurley- 
Keller et al. (1998). The use of modern isochrones that take 
diffusive processes, recent advances in nuclear reaction rates, 
etc., into account can be expected to have some impact on 
age determinations, even if the same distances and chemical 
abundances are assumed. It is perhaps especially encourag¬ 
ing that current stellar models are able to reproduce the mor¬ 
phologies and the locations of the different photometric se¬ 
quences in Carina very well once relatively small adjustments 
are made to the predicted B-V colors. Something worth men¬ 
tioning is that, as shown many years ago by, e.g., Castellani & 
DeglTnnocenti (1995, and references therein), isochrones for 
young ages and very low metallicities have relatively smooth 
turnoffs, in stark contrast with the pronounced blueward hook 
features that are characteristic of their counterparts at near so¬ 
lar abundances. In a qualitative sense, our 2-A Gyr isochrones 
look very similar to those computed > 20 years ago. (Because 
the model grids used in this study are limited to ages ^ 5 Gyr, 
the remainder of this paper will concentrate on just the oldest 
stellar populations of Carina.) 

Figure 6 illustrates the superposition on the Carina CMD 
of 11 and 13 Gyr isochrones for metallicities that span nearly 
the full ranges in [Fe/H] and [a/Fe] that have been derived 
spectroscopically. (Fully consistent ZAHBs, to be discussed 
near the end of this subsection, are also shown.) In confirma¬ 
tion of the remarks made above concerning the oxygen abun¬ 
dance, the isochrones for [Fe/H] = -2.0 and [a/Fe] = 0.0 pre¬ 
dict very similar turnoff and SGB luminosities as those for 
the same ages but for [Fe/H] = -2.40 and [a/Fe] = -1-0.4 (see 
the previous figure). Both assume exactly the same value of 
[O/H] and, as a close inspection of both plots reveals, their 









8 


VandenBerg, Stetson, & Brown 


locations relative to the observed TO and SGB stars in their 
vicinity are nearly the same. However, the RGB portions of 
these isochrones are obviously not coincident, but that can be 
attributed to the differences in the abundances of the heav¬ 
ier metals, notably the iron-peak elements. Varying the oxy¬ 
gen abundance does not affect the location of the RGB at low 
metallicities (see VI2). 

Similarly, considering the isochrones for [Fe/H] = -2.8 and 
[a/Fe] = -tO.4, on the one hand, and those for [Fe/H] = -2.0 
and [a/Fe] = 0.0, on the other, it is readily understood that it is 
the 0.4 dex difference in [O/H] which is primarily responsible 
for the ^ 1 Gyr difference in age at a given TO luminosity (see 
Fig. 6). Although the abundances of the iron-peak elements 
differ by 0.8 dex, and those of the other a-elements (Ne, Mg, 
Si,...) differ by 0.4 dex, these variations mainly cause a small 
shift in the model scale due to their effects on the atmo¬ 
spheric opacities (see V12). The implication of these results is 
that the age at a given TO luminosity varies quite weakly with 
position along the sloped parts of the [a/Fe] vs. [Fe/H] rela¬ 
tion shown in Fig. la. Conversely, a given subgiant thickness 
would imply a narrower range in [Fe/H], at a common age, if 
[a/Fe] = constant than if [a/Fe] and [Fe/H] are anti-correlated 
(as observed). 

This undoubtedly helps to explain why Bono et al. (2010) 
argued that the oldest TO/SGB stars in Carina span a signif¬ 
icantly smaller range in [Fe/H] than that derived from spec¬ 
troscopic work. They reached this conclusion as the result of 
(i) comparing the fiducial sequences for a few of the Galactic 
GCs, which have [a/Fe] r; 0.4, with the Carina CMD, and (ii) 
examining synthetic CMDs that were derived from isochrones 
for [a/Fe] = 0.0 and a wide range in [Fe/H]. Both of these av¬ 
enues assumed, either directly or indirectly, constant [a/Fe]; 
in which case, isochrones of similar age but different [Fe/H] 
(and [O/Fe], since i5[0/H] = (5[Fe/H]) will predict a larger 
difference in their TO/SGB luminosities than one expects if 
[a/Fe] decreases with increasing [Fe/H]. What is of critical 
importance in this discussion is how [O/H] (which is equiv¬ 
alent to the absolute O abundance, if the models being com¬ 
pared assume the same solar abundances) varies with [Fe/H]. 

Perhaps the most compelling argument presented by Bono 
et al. (2010) is that, as demonstrated in their CMD simula¬ 
tions, the component of Carina’s HB that is populated by its 
oldest stars is too thin to be consistent with the range in [Fe/H] 
that has been derived in spectroscopic studies. Indeed, this is 
also suggested by the ZAHB models that appear in the bot¬ 
tom panel of Fig. 5. Most of the old core He-burning stars 
(i.e., those populating the nearly flat portion of the HB at 
(B-V)o ^ 0.55), are brighter than the ZAHB for [Fe/H] = -2.0 
(the short-dashed curve) and fainter than the one for [Fe/H] 
= -2.4 (the solid curve). Moreover, there is no evidence at all 
for old HB stars with [Fe/H] ^ -1.8 as they should lie below 
the ZAHB represented by the long-dashed curve. 

Before moving to the next figure, note that the last of these 
three ZAHBs, which has been extended to a mass that is pre¬ 
dicted for a 2 Gyr RGB tip star, reaches close to the maximum 
luminosity of the HB clump (as shown). Encouragingly, a star 
of the same [Fe/H] that has an RGB tip age of 6 Gyr (approx¬ 
imately the age of the oldest lA stars) is predicted to reach 
its ZAHB location at close to the minimum luminosity of the 
clump (at the location of the long-dashed curve). Accord¬ 
ing to our models, and consistent with the observed MDF to 
within its uncertainties, few of the observed HB clump stars 
are predicted to have [Fe/H] >-1.4 (which has been assumed, 
together with [a/Fe] = 0.0, in the computation of the ZAHB to 


which a small filled circle has been attached at its bright end), 
since this ZAHB matches the red edge of the distribution of 
HB clump stars. (These abundances are consistent with the 
variation of [a/Fe] with [Fe/H] given by the solid curve in 
Fig. la.) To reinforce this point, a ZAHB for [Fe/H] = -1.2 
and [a/Fe] = 0.0 would have been located well to the red of 
the HB clump, almost intersecting the observed RGB, had we 
chosen to plot this case. (There may well be some stars in 
our CMD that have [Fe/H] ^ -1.2, in accordance with spec¬ 
troscopic findings, but the small minority of such stars pre¬ 
sumably does little more than contribute to the scatter in the 
vicinity of the HB.) 

The benefits of being able to compare observations with 
stellar models for the observed chemical abundances is fur¬ 
ther revealed by the ZAHBs that have been plotted in Fig. 6. 
They show that the small luminosity width of Carina’s old 
HB does not preclude the presence of stars with [Fe/H] val¬ 
ues that vary by more than 1 dex, provided that [a/Fe] varies 
with [Fe/H] in approximately the observed way (i.e., follow¬ 
ing the trend given by the dashed line in Fig. la). The ZAHB 
for [Fe/H] =-1.4 and [a/Fe] = -0.4 provides quite a good fit 
to the faintest of the cluster HB stars that are fainter than the 
clump stars and that have 0.35 ^ (B-V)o ^ 0.55. This is not 
unexpected since, as reported by VandenBerg et al. (2000), 
the luminosity of a ZAHB is quite a strong, inverse function 
of [a/Fe], at fixed values of [Fe/H] and TQ 

To be sure, the distribution of the stars relative to the three 
ZAHBs that have been plotted in Fig. 6 seems anomalous. 
Most of the observed HB stars are fainter and redder than the 
ZAHB for [Fe/H] = -2.0 and [a/Fe] = 0.0, implying that they 
are more metal rich, which is obviously in conflict with the 
MDF. A possible solution to this dilemma is that the oxygen 
abundance is significantly higher than [O/Fe] = -tO.4 in (at 
least) the most metal-deficient stars in Carina, and that it re¬ 
mains offset from the [m/Fe] values of the other a-elements 
by ^ 0.2 dex (or more?) as [a/Fe] decreases with increas¬ 
ing [Fe/H]. For instance, [O/Fe] ^ -tO.6 is generally found 
in Milky Way halo stars (e.g., Fabbian et al. 2009, Ramirez 
et al. 2012, V14b), whereas a smaller enhancement ([m/Fe] 
= 0.3-0.4) is generally found for Mg, Si, Ca, Ti, etc. (see, e.g., 
Cayrel et al. 2004). At low metallicities, a ZAHB is stretched 
to redder colors and fainter luminosities over a considerable 
fraction of its length as the absolute oxygen abundance is in¬ 
creased (see, e.g., VandenBerg & Bell 2001; their Fig. 6). (At 
higher [Fe/H], enhanced [O/H] mainly causes a redward dis¬ 
placement of just the red end of a ZAHB.) Thus, higher oxy¬ 
gen could shift the ZAHB loci for [Fe/H] ^ -2.0 in just the 
way that is needed to reconcile them with the locations of a 
signficant fraction of the reddest stars in the “old" HB compo¬ 
nent of Carina. Unfortunately, little seems to be known about 
the [O/H] values in Carina stars, and it must be left for future 

® Work is in progress by D.A.V. to produce ZAHB loci for all of the chem- 
ical abundance choices considered by V14a, and we leave a detailed dis¬ 
cussion of those ZAHBs to the forthcoming paper. However, it may be of 
some interest to note that, despite the revisions which have been made to 
the Victoria evolutionary code over the past 15 years, current ZAHB mod¬ 
els predict very similar luminosities for RR Lyrae stars as those presented 
by VandenBerg et al. (2000). Whereas the latter neglected diffusion physics 
and assumed a primordial helium abundance corresponding to 7 = 0.235, 
which was the favored value at that time, current Victoria models take the 
gravitational settling of helium into account while assuming the current best 
estimate of the primordial value of Y 0.250; e.g., Komatsu et al. 2011). 
Fortuitously, this increase in Y largely compensates for the effects of diffusion 
on the properties of RGB tip stars (primarily the envelope helium abundance) 
that control the predicted luminosity of the HB. 
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Figure 7. Best-fit model (grey shading) for the oldest Carina stars (blue 
points), assuming a single burst of star formation and the [Fe/H] distribution 
of Figure lb. The relative weights (VF) are also specified in the inset table 
for each panel. The o-varying model produces a much narrower SGB than 
observed {left panel). The O-varying model (with [a/Fe] otherwise held at 
-1-0.4) produces a somewhat broader SGB that is still naiTower than observed 
(right panel). 

work, both observational and theoretical, to investigate these 
speculations. 

Figure 6 also suggests that the faintest SGB stars are old, 
iron-rich, a-poor stars. (The delay time between SF and the 
onset of SN Type la explosions is not necessarily an argument 
against the near coevality of stars with very different [a/Fe] 
values, as it can be as short as ^ 100-500 Myr or as long as a 
few Gyr, depending on, among other things, whether the star 
formation rate is high or low, respectively; see the review by 
Maoz et al. 2014.) It is possible that some of them are field 
stars, though we have tried to remove such objects. Alterna¬ 
tively, it may be that the actual [O/Fe] values are higher than 
those implied by the assumed values of [a/Fe]. This would 
have the effect of reducing the predicted age at a given TO lu¬ 
minosity, or equivalently, of shifting isochrones for a fixed age 
and [Fe/H] to higher My. That is, the faintest SGB stars could 
be very old stars with low iron, and relatively high oxygen, 
abundances (that evolve to red ZAHB locations). To provide 
a limited investigation of the consequences of this possibility 
for the oldest TO/SGB stars in Carina, some of the simula¬ 
tions presented in the next subsection will make use of the 
model grids for [a/Fe] = -1-0.4 and high values of [O/Fe] that 
were employed in the recent study of several UFD galaxies by 
Brown et al. (2014). 

4.2. Simulations of the SFH of the Oldest Stellar 
Populations 

To explore the SFHs that are consistent with the oldest pop¬ 
ulations in Carina under different metallicity assumptions, 
we have constructed several synthetic CMDs in the B and 
V bands, employing the same calibrated isochrone library 
described above. Two different SFH scenarios have been 
explored: a short burst of star formation with an idealized 
MDF (Fig. lb), and an extended period of star formation last¬ 
ing several Gyr, with only loose constraints on the MDF. In 
each of these two SFH scenarios, the relationships for [a/Fe] 
vs. [Fe/H] and [O/Fe] vs. [Fe/H] were also varied. A linear 
combination of synthetic CMDs, each representing a partic¬ 


ular age and metallicity, was used to fit the observed CMD. 
These fits were performed via the minimization of the Pois¬ 
son maximum likelihood statistic of Dolphin (2002), and the 
quality of the best fit was evaluated by comparing the best- 
fit score to the distribution of such scores arising from Monte 
Carlo realizations of the photometric data. 

Usually, synthetic CMDs are constructed using the results 
of artificial star tests. Such tests characterize the photomet¬ 
ric uncertainties and incompleteness as a function of position 
within the CMD, allowing each point along an isochrone to 
be scattered in the synthetic color-magnitude plane. Due to 
the heterogeneous nature of the observations studied here, ar¬ 
tificial star tests are impractical, but the photometry of each 
star in the catalog comes from many individual exposures, 
such that the photometric uncertainties are well understood. 
Furthermore, for the CMD region being simulated, the photo¬ 
metric errors are small (<0.05 mag); consequently, the com¬ 
pleteness should be near 100%. For these reasons, we used 
the photometric errors from the catalog to define a Gaussian 
scattering kernel as a function of position within the CMD. 

The synthetic CMDs were populated using a Salpeter 
(1955) initial mass function (IMF), although the small mass 
range of stars in the fitting region (< 0.1 Ad©) makes the as¬ 
sumed IMF unimportant. We assumed a binary fraction of 
14% (Minor 2013), but the results would not change signif¬ 
icantly if we were to adopt a value closer to 50%, which is 
typically found in other dwarf galaxies. To fit the SFH of the 
oldest population, we isolated the faintest SGB stars in the 
catalog, along with its neighboring stars in the vicinity of the 
MSTO and the lower RGB (see Figure 7). 

We begin with a simple model in which (i) the star forma¬ 
tion is assumed to occur in a single burst, and (ii) the variation 
in the number of stars as a function of [Fe/H] corresponds to 
the MDF shown in Fig. lb. However, we have explored two 
distinct assumptions regarding [a/Fe] and [O/Fe]. The first is 
the “a-varying” model, where [a/Fe] increases with decreas¬ 
ing [Fe/H] (as in Figure la). The second is the “O-varying” 
model, where [O/Fe] increases with decreasing [Fe/H], while 
all other a elements are held at a constant [a/Fe] = h- 0.4 (as 
typically found in old, very metal-poor populations). The 
adopted relation between [O/Fe] and [Fe/H] is the same as 
that assumed by Brown et al. (2014), which is based on data 
presented by Frebel (2010). 

Figure 7a shows the best-fit burst with the a-varying model. 
The simulated SGB is significantly narrower than the ob¬ 
served SGB, resulting in a rather poor fit (Xeff= 1 ■8)- Figure 7b 
shows the best-fit burst in the case of the O-varying model. 
The simulated SGB is broader, but still not as broad as that 
observed, again producing an unsatisfactory fit (Xgg=1.4). Al¬ 
though not shown here, we note that a model in which all of 
the a-elements, including oxygen, have a constant enhance¬ 
ment (specifically, [a/Fe] = -1-0.4) for all [Fe/H] values looks 
quite similar to that illustrated in Figure 7b. Both of these fits 
are ruled out at more than 6 cr. Thus, if the MDF in the oldest 
population of Carina is similar to that shown in Figure lb, it 
is unlikely to be a single-age population. 

We now turn to a model with a more complex SFH, shown 
in Figure 8. Instead of a single burst, the age distribution is 
only restricted to the range of 8 to 14.4 Gyr. The MDF is 
constrained to iron abundances -3.0 < [Fe/H] < -1.4, as in 
the case of the single burst models, but it is permitted to take 
any form within this range. (These fits were performed using 
the StarFISH code of Harris & Zaritsky 2001.) Our results 
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Figure 9. As in Figure 8, but for an assumed distance modulus that is 0.1 mag larger. 
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for an a-varying model and an O-varying model are shown 
in the upper and lower panels, respectively. For both cases, 
Xeff=l-2 (rounded to one decimal place), but the O-varying fit 
is the slightly better one. However, although nearly the same 
value of the goodness-of-fit statistic is obtained, the derived 
weighting of the isochrones in Fig. 8d resembles the observed 
MDF more so than that shown in Fig. 8c. In particular, the 
former simulation predicts more stars with [Fe/H] ^ -2.2 (re¬ 
call Fig. lb) than in lower or higher metallicity bins, whereas 
the opposite is found in the latter simulation. 

It may be recalled that the results presented in Fig. 8c were 
anticipated by the overlays of isochrones to the Carina CMD 
discussed previously. As shown in Fig. 6, isochrones for high 
ages and [Fe/H] values are apparently needed to explain the 
faintest SGB stars if the assumed a-element abundances fol¬ 
low the adopted relation between [a/Fe] and [Fe/H]. As indi¬ 
cated by the simulation results in Fig. 8d, this is not the case 
if the stellar models allow for high oxygen abundances. Also 
worth mentioning is the fact that significantly higher ages, at 
a fixed metallicity, would have been found had we neglected 
diffusive processes (see, e.g., VandenBerg et al. 2002, their 
Fig. 2). (However, if such physics is ignored, it would be¬ 
come especially difficult to reconcile stellar ages with the age 
of the universe; see V14b.) 

Although the fits are not perfect, both of the cases con¬ 
sidered in Fig. 8 demonstrate that the oldest population in 
Carina can be reproduced if the star formation epoch lasted 
a few Gyr. Although there is a general trend of decreasing 
age at increasing metallicity, the particular distribution of the 
points in the age-metallicity plane depends quite sensitively 
on how [O/Fe] varies with [Fe/H], and the distance modulus 
assumed. A shorter distance would be more problematic than 
a longer distance because it would imply that all of the stars 
(notably those on the SGB) have fainter My values, thereby 
making them older. As a consequence, the weighting of the 
isochrones found for the O-varying model (Fig. 8d) would 
move in the direction of that obtained for the a-varying model 
(Fig. 8c), and the results for the latter case would be skewed to 
even older ages and higher [Fe/H] values (in greater conflict 
with the observed MDF). 

Figure 9 shows that the opposite trend is found if a slightly 
larger distance modulus (by 0.1 mag) is assumed. Compar¬ 
ing, e.g.. Fig. 9c with Fig. 8c reveals that the distribution and 
weighting of the points has been shifted to lower [Fe/H] val¬ 
ues. That is, the models predict larger/racf/oni of older, more 
metal-poor stars if the distance to Carina is 0.1 mag larger 
than our adopted (ZAHB-based) modulus. (The same thing, 
which may seem counter-intuitive, is implied by the differ¬ 
ences between Figs. 9d and 8d.) As a consequence, the a- 
varying model provides a more agreeable fit to the observa¬ 
tions, while the O-varying model seems to be somewhat less 
satisfactory. However, the main value of the simulations pre¬ 
sented in Figs. 8 and 9 is to illustrate how they are affected 
by parameter variations. To obtain accurate and robust age 
determinations, it is clearly of paramount importance to have 
reliable abundances of the individual a-elements (especially 
O, but also Mg and Si; see V12) so that their variations with 
[Fe/H], as well as the observed star-to-star scatter of [m/Fe] at 
a fixed iron abundance, can be taken into account when simu¬ 
lating the observed CMD of Carina. 

4.3. The Thin RGB of Carina 

The fact that Carina has a very narrow RGB is unexpected 
given the wide ranges in ages and metal abundances encom- 



0.4 0.8 1.2 0.8 1.2 1.6 


(B-V)„ (V-I)„ 

Figure 10. Fits of the RGB segments of 13.0 Gyr, [Fe/H] = -2.40 and 12.5 
Gyr, [Fe/H] = —1.55 isochrones to the giant branches of M92 and M3, re¬ 
spectively. (M92 has a bluer RGB than M3.) Both isochrones assume 
Y = 0.25 and [ci/Fe] = 0.40 and both were adjusted by S{B-V) = -0.01 and 
<5(y — 7) = -0.03 mag in order to provide the best overall fits to the cluster 
observations. 

passed by its stellar populations. As Venn et al. (2012) have 
remarked in their recent study of this dSph, it is possible that 
the giants “have a fortuitous alignment in the age-metallicity 
degeneracy, such that the older metal-poor stars overlie the 
metal-enhanced intermediate-age stars". Although this over¬ 
lap is undoubtedly happening, to what extent does the ob¬ 
served RGB constrain the properties of Carina’s stellar pop¬ 
ulations? In what follows, we will endeavor to answer this 
question using several of the same isochrones that seem to 
provide quite good fits to the TO and HB observations. 

Fortunately, the predicted morphologies of current Victoria- 
Regina isochrones reproduce the shapes of GC giant branches 
very well when standard reddenings (from Schlafly & 
Finkbeiner 2011), distances based on ZAHB models that sat¬ 
isfy RR Lyrae constraints (see V13, V14b), and metallici- 
ties close to those derived by, e.g., Kraft & Ivans (2003) and 
Carretta et al. (2009a) are assumed. As shown in the left- 
hand panel of Figure 10, the isochrones for the same ages 
and chemical abundances that were fitted to the turnoff obser¬ 
vations of M 3 and M 92 in § 13.11 on the assumption of the 
cluster properties that are given explicitly in Fig. 3, provide a 
close match to their RGB populations on the [(B-y)o, My]- 
diagram. To accomplish this, the B-V colors from the CV14 
transformations had to be corrected by just -0.01 mag, inde¬ 
pendently of gravity and T^ff. The right-hand panel shows that 
equally good agreement can be obtained on the [(y-/)o, My]- 
diagram; but in this case, the required color adjustment turns 
out to be S(V-I) = -0.03 mag. (The origin of the small zero- 
point offsets is not known, but they can be reasonably at¬ 
tributed to small photometric errors or minor problems with, 
e.g., the CV14 color transformations, the model Teff scale, 
and/or the assumed cluster properties.) 

What is especially encouraging is that exactly the same 
6(B - V) and 6(V - /) color offsets are needed for both M 3 
and M 92, despite the difference of nearly 1 dex in their [Fe/H] 
values. This implies that models which are fitted to the giant 
branch of Carina should be trustworthy in an absolute sense 
if we simply apply the same color adjustments to them. In 
effect, we have calibrated the RGBs of the Victoria-Regina 
isochrones using M92 and M3. (As shown by V14a, these 
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Figure 11. Left-hand panel: Overlay of the RGB portions of three isochrones onto the Carina CMD. The bluer of the two solid curves assumes [Fe/H] = -2.40, 
[o/Fe] = 0.40, and an age of 13 Gyr; the other one was generated for [Fe/H] = —1.30, [a/Fe] = 0.07, and an age of 6 Gyr. The dashed curve is the best-fit isochrone 
to the M3 RGB (see the previous figure). Middle panel: As in the left-hand panel, except that the dashed curve has the indicated properties and the thin solid 
curves represent the RGB portions of 4, 5, and 6 Gyr isochrones for [Fe/H] = -1.8 and [a/Fe] = 0.40. Right-hand panel: The same 6 and 13 Gyr isochrones that 
appear in the left-hand and middle panels have been replotted, along with three isochrones for the indicated ages and chemical properties (the thin solid curves). 
The adopted reddening and distance modulus of Carina are the same as in previous figures, and we have assumed that EiV—1) = \.2AAE{B — V) (Casgrande & 
VandenBerg 2014). The predicted colors along the isochrones were adjusted by 5(y—I) = -0.03 mag prior to plotting them. 


computations are equally successful over a much broader 
range in [Fe/H]. The main difficulty, as also reported by V13, 
is that somewhat different offsets are needed to fit the MSTO 
observations than to match the location of cluster RGBs; i.e., 
the predicted difference in color between the TO and the lower 
giant branch differs from the observed difference by a small 
amount. However, this is not of particular concern because, as 
noted at the end of the previous paragraph, such discrepancies 
could easily be due to relatively minor deficiencies in some of 
the more uncertain aspects of the stellar models, the assumed 
chemical compositions, etc.) 

It is already apparent in Fig. 10 that a difference of 0.85 dex 
in [Fe/H] implies a much thinner RGB on the [(y-/)o, My]- 
diagram, especially at My ^ 0.0, than on the [(B-y)o, My]- 
plane. This is the expected consequence of the fact that the 
B-V color index has a much stronger metallicity dependence 
than V -I. Because the latter will also be less affected by 
variations in the metals mixture, it should be more straightfor¬ 
ward to interpret the fits of isochrones to the y/ photometry of 
Carina than to CMDs based on BV data. Therefore, we will 
hrst examine the implications for the observed RGB on the 
[(y-/)o, My]-diagram of the isochrones that have been com¬ 
pared with the TO observations in the previous section. {V-I 
colors were not considered previously because the brightness 
of the night sky in I makes it impossible to obtain ground- 
based photometry in this filter for the faintest TO stars in Ca¬ 
rina.) 

Figure 11 compares the giant branches of several 
isochrones with our VI data, assuming that Carina has E(B- 
y) = 0.054 and {m—M)y = 20.28, as adopted throughout this 
investigation. To deredden the V-I colors, we have assumed 
that E{V-I) = 1.244£(B-y) (see CV14). Since > 95% of 
the stars in the dSph appear to have [Fe/H] values between 
-2.8 and -1.3 (according to the MDF in Fig. lb), we have 


plotted (in all three panels, as thick solid curves) the RGBs of 
a 13 Gyr isochrone for [Fe/H] = -2.8, [a/Fe] = +0.4 and a 6 
Gyr isochrone for [Fe/H] =-1.3, [a/Fe] =+0.07 (consistent 
with the relation between [a/Fe] and [Fe/H] shown in Fig. la). 
It is apparent (see the left-hand panel) that these isochrones 
bracket the upper-RGB populations quite well and that they 
predict an exceedingly thin lower RGB. The dashed curve 
represents the giant-branch segment of the the same isochrone 
that was superimposed on the CMD of M 3 in Fig. 10. This 
is signihcantly redder than the majority of the bright giants in 
Carina, from which one would be inclined to conclude, in the 
absence of any additional information, that the galaxy has few, 
if any, stars with [Fe/H] values as high as -1.55. However, we 
know that from spectroscopic work that there are large num¬ 
bers of giants with higher iron, but lower a-element, abun¬ 
dances, and photometric studies have revealed the existence 
of many young stars. In accordance with standard stellar evo¬ 
lutionary theory (e.g., see V14a), the most iron-rich giants 
in Carina are bluer than the RGB of M 3, at the same My, 
because the former are younger and a-poor (in an absolute 
sense) relative to their cluster counterparts Q 

In the middle panel of Fig. 11, a 12 Gyr isochrone for 
[Fe/H] = -1.6 and [a/Fe] = -0.27 has been overlaid onto 
the Carina CMD, along with the same two isochrones from 
the left-hand panel (the thick solid curves) and 4-6 Gyr 
isochrones for [Fe/H] =-1.8, [a/Fe] =+0.4 (the thin solid 
curves). The separation in color at a fixed magnitude between 
the old -2.8 and -1.6 isochrones is quite small — easily con¬ 
sistent with the observed spread in RGB colors at a given lu¬ 
minosity (at least on this CMD). Whereas the three isochrones 
for intermediate ages lie between those for [Fe/H] = -2.8 and 

^ It should be kept in mind that it is the absolute abundances of the metals, 
not [m/Fe] ratios, that are important for stellai' models. 
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Figure 12. As in the previous figure, except that the comparisons between theory and observations are made on the [(B — V)q, Mv]-diagYam and the isochrone 
colors were corrected by 5(B — V) = -0.01 mag. 


-1.3 at My ^ -1, they are located slightly to the blue of the 
lowest metallicity isochrone at fainter absolute magnitudes. 
Even so, the width of the RGB below the HB clump is still 
very narrow, and it would remain so even if the range in 
[Fe/H] or age were expanded somewhat. In fact, the predicted 
RGB is actually a bit thinner than the observed one, which 
may be a consequence of photometric errors, small as they 
are (see §|2]l. Note that old (~ 12 Gyr) isochrones for [Fe/H] 
< -1.6 would overlie the 4-6 Gyr isochrones that have been 
plotted for [Fe/H] =-1.8, which illustrates the well-known 
age-metallicity degeneracy. 

Very similar comparisons between theory and observa¬ 
tions are shown in the right-hand panel; in this case, the 
intermediate-age isochrones are for [Fe/H] =-1.4 and [a/Fe] 
= -1-0.13. Not unexpectedly, they lie just to the blue of the 6 
Gyr isochrone for [Fe/H] = -1.3 and [a/Fe] = -1-0.07. A sim¬ 
ilar set of 4-6 isochrones was computed for [Fe/H] = -1.6 
and [a/Fe] = -1-0.27, which is also consistent with the [a/Fe] 
vs. [Fe/H] relation in Fig. la, and although we have chosen 
not to plot them here, they lie approximately midway between 
those for [Fe/H] = -1.8 (middle panel) and -1.4 (right-hand 
panel). As a result, they would overlie the densest concentra¬ 
tion of Carina giants better than the latter, which is perhaps 
to be expected since the peak in the MDF occurs at [Fe/H] 
~ -1.65 (see Fig. lb). 

In hindsight, it is not too much of a surprise that the Ca¬ 
rina RGB is very narrow on the [(y-/)(), My]-diagram, given 
that the vast majority of its stars appear to be more metal-rich 
than M 92, more metal-poor than M 3 (largely due to the low 
abundances of the a-elements), and older than ~ 3 Gyr (see 
Fig. 5). Because nearly all of the giants in the dSph lie to the 
blue of M 3’s RGB, (see the left-hand panel of Fig. 11), it is to 
be expected that the thickness of the galactic RGB will be less 
than the separation between the M 92 and M 3 giant branches, 
which amounts to only 6{V-I) Ki 0.04 mag at My values that 
are fainter than the cluster HBs (see Fig. 10). Young (< 3 
Gyr), [Fe/H] ^-1.8 Gyr giants would be located on the blue 
side of a 13 Gyr, [Fe/H] = -2.8 isochrone, but they do not ap¬ 


pear to be present, at least in significant numbers (see Fig. 11). 
We conclude that the thin RGB of Carina is completely con¬ 
sistent with the observed and derived properties of its stellar 
populations, even if moderately large uncertainties are taken 
into account. 

However, V-I colors are a better discriminator of differ¬ 
ences in Teff than of metallicity, so it is instructive to perform 
the same comparisons of isochrones with the RGB of Carina 
as those just discussed, but on the [(B-Vjo, My]-diagram. In¬ 
deed, the analogous plot (see Figure 12) differs by more than 
one might have anticipated from the previous one, though 
there are several similarities as well. The most noticeable dif¬ 
ference is that the giants are not as centrally located between 
the two thick solid curves as in Fig. 11 (compare the left-hand 
panels). Indeed, the 6 Gyr isochrone for [Fe/H] = -1.3 and 
[a/Fe] = -1-0.07 is significantly redder than the red edge of the 
distribution of giants, though a much better fit to the latter is 
obtained by a 4 Gyr isochrone for the same metal abundances 
(judging from the location of the isochrone for [Fe/H] =-1.4 
and [a/Fe] = -1-0.13 in the right-hand panel). In addition, the 
median RGB of Carina seems to be matched quite well by 
5-6 Gyr isochrones for [Fe/H] = -1.8 and [a/Fe] = -1-0.4 (see 
the middle panel). Both of these isochrones were located just 
to the left of the main concentration of bright giants on the 
[(y-/)o,Mv]-plane. 

As found previously, a 13 Gyr isochrone for [Fe/H] = -2.8 
and [a/Fe] = -1-0.4 provides a good fit to the blue edge of the 
distribution of Carina giants. Why the same interpretation 
of the redder stars is not found from both B-V and V-I 
colors is not clear, but there are a few possible explanations 
for this. For one, it is possible that the [Fe/H] scale derived 
by Starkenburg et al. (2010), which has been assumed here 
(see Fig. lb), is too metal-rich by ^ 0.15-0.2 dex. Such a 
shift, which is < the differences in measured GC metallici- 
ties over the years (e.g., see Carretta & Gratton 1997, Kraft 
& Ivans 2003, Carretta et al. 2009a) would affect the hori¬ 
zontal separation between the two isochrones in the left-hand 
panel of Fig. 12 considerably more than in the case of the 
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corresponding isochrones in Fig. 110 We see, for instance, 
that a 12 Gyr isochrone for [Fe/H] = -1.8 provides a much 
better fit to the red edge of the main distribution of Carina gi¬ 
ants (right-hand panel) than the 12 Gyr isochrone for [Fe/H] 
= -1.6 that has been plotted in the middle panel. (If it is only 
the iron abundance that changes, the derived value of [a/Fe] 
would increase by the amount that the value of [Fe/H] is re¬ 
duced. However, one could compensate for that by decreas¬ 
ing the assumed age. For example, a 10 Gyr isochrone for 
[Fe/H] = -1.8 and [a/Fe] = -1-0.02 is essentially indistinguish¬ 
able from one for the same [Fe/H], but [a/Fe] =-0.13; see the 
right hand panel.) 

This brings us to the second possibility, which is that stars 
of different metal abundances have different ages (as sug¬ 
gested by our simulations in the case of the oldest stars). 
Fig. 12 shows that a 6 Gyr isochrone for [Fe/H] = -1.8 and 
[a/Fe] = -1-0.4 (middle panel) is bluer than a 4 Gyr isochrone 
for [Fe/H] = -1.4 and [a/Fe] = -1-0.13 (right-hand panel) by 
only a small amount. In fact, the latter would be nearly coin¬ 
cident with the former if it had an age of 3 Gyr. (This expla¬ 
nation would indeed imply that ages vary with metallicity in 
just the right way to produce a thin giant branch.) In addition, 
variations in the heavy-element mixture may explain some of 
the differences in the fits of isochrones to the BV amd VI pho¬ 
tometry of Carina giants. It is an approximation to adopt the 
same abundances for all of the a-elements at a given [Fe/H] 
value. As shown by VI2, the location of the RGB on the 
theoretical plane has a comparable dependence on the abun¬ 
dance of Si as of Mg, and if, e.g., [Si/Fe] < [Mg/Fe], then the 
values of [a/Fe] that we have assumed are too high, and the 
isochrones (especially at higher metallicities) are too red. Un¬ 
fortunately, not a great deal is known about the abundances of 
silicon in Carina. (Star-to-star abundance variations will also 
have important consequences for stellar temperatures and col¬ 
ors, but color transformations that allow for differences in the 
abundances of single elements are not yet available.) At this 
stage, the uncertainties are still large enough that the minor 
inconsistencies that exist between Figs. 11 and 12 are not of 
particular concern. 

5. SUMMARY 

The most important result of this investigation is that stellar 
models are able to reproduce the main features of the CMD 
of Carina quite well when they assume the full range in metal 
abundances (^1.5 dex in [Fe/H]) that has been derived from 
high-quality spectroscopic data. This contradicts previous 
findings, notably by Bono et al. (2010), who argued from their 
analysis of the observed photometry that the metallicity vari¬ 
ations encompassed by either the oldest or intermediate-age 
populations cannot amount to much more than a few tenths of 
a dex. The key to our success is that we have used isochrones 
(from V14a) which assume heavy-element abundances that 
follow the observed relations between [a/Fe] and [Fe/H]. (In¬ 
terpolations in those grids can be made for any value of [a/Fe] 
between -0.4 and -1-0.4 at a fixed [Fe/H] value.) In fact, it 
would not have been possible to reproduce some CMD fea¬ 
tures, such as the small range in luminosity spanned by the 

* If anything, spectroscopic studies in the past few years have caused even 
more consternation regarding absolute [Fe/H] determinations than earlier 
ones. For instance, whereas Carretta et al. (2009a) found that M 15, M92, 
and NGC5466 all have [Fe/H] values within the range of -2.31 to —2.35, 
Sobeck et al. (2011) and Roederer & Sneden (2011) have obtained [Fe/H] 
;$ —2.6 for M 15 and M92, respectively, while Lamb et al. (2015) have de¬ 
rived a value of-1.97 for NGC 5466. 


old HB component of Carina, using models for the measured 
range in [Fe/H] unless [a/Fe] is a strongly decreasing function 
of the iron abundance (as observed). In other words, models 
for constant [a/Fe], irrespective of the value that is adopted, 
are incapable of providing a satisfactory explanation of the 
Carina CMD. 

Being able to transform our models to the [{B-V)o,Mv]- 
and [(y-/)o, Mi/]-diagrams using fully consistent transfor¬ 
mations (from CV14) has certainly made a difference as well, 
because predicted BVI magnitudes have a considerable de¬ 
pendence on the abundances of the a-elements at fixed values 
of [Fe/H], logg, and Teff. Furthermore, we have taken the 
additional step of determining, and applying, corrections to 
the CV14 (B-y)-reff relations in order that isochrones repro¬ 
duce the CMDs of M 3 and M 92, when very good estimates 
of their reddenings, metallicities, and distances are assumed. 
This “calibration" of the isochrones ensures that the models 
faithfully reproduce the locations of the MSTOs and lower 
RGBs, as well as the slopes of the SGBs in GCs that span a 
range of nearly 1 dex in [Fe/H], thereby giving us confidence 
that they will accurately represent the properties of stars in 
the dSph. This procedure should minimize errors that would 
otherwise be present, and possibly be quite substantial. 

The veracity of this assertion is easily verified. There is a 
general tendency of Victoria-Regina isochrones to suffer from 
small zero-point and systematic offsets when compared with 
GC CMDs, though the predicted slopes of the MS and RGB 
look fine (see, e.g., the plots provided by V13 and V14a). 
When compared with the same HST observations that were 
the subject of the paper by VI3, Dartmouth isochrones appear 
to be considerably more problematic (see Dotter et al. 2010; 
the right-hand panels of their Figs. 4 and 5). Differences in 
the color transformations that were employed in the respective 
studies are likely to be mostly responsible for this, in view of 
the fact that very similar physics has been incorporated in both 
the Victoria and Dartmouth codes. However, the same cannot 
be said of BASTI isochrones (e.g., Pietrinferni et al. 2004) 
because these computations do not take diffusive processes 
into account, which have significant consequences for turnoff 
luminosity versus age relations and the predicted Tgff scale, 
especially at low metallicities and high ages (see VandenBerg 
et al. 2002; their Fig. 2). 

It is risky to use stellar models from any source with¬ 
out carefully examining how well they satisfy empirical con¬ 
straints. This is especially true when attempting to interpret 
the photometry of a galaxy, like Carina. When an isochrone 
is overlaid onto the observed CMD, it should superimpose 
stars of nearly the same age and chemical abundances whether 
they are located near the TO, along the SGB, or on the giant 
branch. This is, of course, complicated by known degenera¬ 
cies between the effects of age and chemical abundances, but 
if the adopted isochrones do not reproduce the locations and 
morphologies of GC CMDs particularly well, those models 
are unlikely to provide very meaningful results when applied 
to galactic observations. (GC ages, in contrast, mainly involve 
fits of isochrones to turnoff stars.) 

The ages of the different stellar populations in Carina that 
we obtained from overlays of isochrones onto the observed 
CMD are quite similar to those found previously, especially 
in a relative sense. The oldest of the intermediate-age stars 
appear to have formed 6-6.5 Gyr ago, whereas the range 
in age needed to reconcile stellar models with the faintest 
turnoff/SGB stars is ^ 9-14 Gyr. In agreement with earlier 
determinations, there was apparently a ^ 3-4- Gyr hiatus in 
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the formation of stars in Carina after the first SF epoch. 

Our CMD simulations, which are based on isochrones that 
have been carefully calibrated, rule out the possibility that 
the faintest TO/SGB stars in Carina are a single-age popu¬ 
lation by more than 6 a. However, simulations in which ages 
are limited to the range 8-14.4 Gyr and iron abundances to 
-3.0 < [Fe/H] < “1-4, without constraining how the number 
of stars varies with metallicity, provide a significantly bet¬ 
ter match to the observed CMD if the earliest star forma¬ 
tion epoch lasted a few Gyr. Moreover, judging from the 
derived weighting of the points in a grid of age and [Fe/H], 
there are two ways of obtaining a reasonable approximation 
to the observed MDF (both of which suggest that ages de¬ 
crease somewhat with increasing [Fe/H]). Either the adopted 
distance modulus, (m-M)v = 20.28, should be increased by 
^0.1 mag, if [a/Fe] varies with [Fe/H] according to the re¬ 
lation shown in Fig. la, or the assumed O abundances should 
be increased and the dependence of [O/H] on [Fe/H] is similar 
to that adopted by Brown et al. (2014). In support of the lat¬ 
ter possibility, we note that [O/Fe] ^ -1-0.6 is typically found 
in very metal poor. Galactic halo stars (e.g., Garcia P erez et 
al. 2006, Ramirez et al. 2012), and as discussed in § 14.11 it 
may be difficult to reconcile the distribution of stars along the 
old HB component of Carina with the observed MDF unless 
the stars have relatively high O abundances. (A larger distance 
seems unlikely in view of the recent empirical determinations 
by Pietrzyhski et al. 2009 and Coppola et al. 2013.) 
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Although our isochrones had no difficulty reproducing the 
narrow RGB of Carina on the V -I, V CMD, the observed 
width on the B-V,V color plane was found to be smaller 
than expected. Improved consistency would be obtained if the 
age-metallicity relations that describe the stellar populations 
in Carina have significant slopes, in the expected sense that 
more metal-rich stars are younger. In fact, this is suggested 
by the simulations that we, and de Boer et al. (2014b), have 
carried out for the oldest, and intermediate-age, populations 
of Carina, respectively. (Our simulations were limited to the 
oldest stars because extensions of the V14a grids to younger 
ages are not available.) In addition, or alternatively, this dif¬ 
ficulty may be a reflection of current uncertainties associated 
with the observed [Fe/H] scale, the assumed [a/Fe] values 
(which are based entirely on Mg, as little is known about the 
abundances of Si), and/or the color transformations. While 
a resolution of this problem must be left for future work, it 
is very encouraging that modern stellar models, on the as¬ 
sumption of the chemical abundances which have been de¬ 
rived from high-resolution spectroscopy, are able to explain 
most of the features of Carina’s CMD quite well. 
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